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摘要 :利用 大 气 和 水 文 模型 定量 描述 陆 表 相关 变量 变化 规律 一 直 是 大 气 科 学 和 水 文学 界 的 研究 热点 。 狼 而 ,由 于 我 国 西部 地 区 
站 点 匮乏 ,传统 气象 观测 站 点 已 不 能 满足 大 尺度 地 表 分 量 高 精度 模拟 分 析 的 需求 。 建 立 SWAT 模型 中 国 大 气 同 化 驱动 数据 集 
(China Meteorological Assimilation Driving Datasets for the SWAT model, CMADS) 驱动 SWAT(Soiláànd Water Assessment Tool) fi 
型 (简称 为 CMADS+SWAT 模式 ) ,选取 传统 气象 站 点 稀缺 的 新 疆 精 博 河流 域 为 并 区 , 兄 成 流域 各 地 表 分 量 (如 土壤 湿度 . 雪 深 、 
融雪 ) 校准 .验证 及 其 时 空 关 系 提取 与 分 析 。 分 析 发 现 :CMADS 数据 集 可 很 好 地 驱动 KE SWAT 模式 完成 本 地 化 工作 。 其 
中 ,CMADS+SWAT 模式 在 月 尺度 上 总 体 NSE 效率 系数 均 在 0.659—0.942, 目 尺 度 也 均 在 0.526 一 0.815。 对 流域 内 土壤 湿度 和 
融雪 过 程 进 行 相关 分 析 发 现 : 精 博 河流 域 土壤 湿度 在 年 内 3 一 4 月 份 达到 其 第 宇 次 峰值 ,主要 贡献 来 自 于 流域 内 高 山 融 雪 现 象 ; 
融雪 期 结束 后 ,流域 降水 量 增 加 ,伴随 气温 上 升 等 现象 导致 土壤 温度 旺 现 波动 态势 ,至 10 月 中 旬 冷 空气 过 境 产 生 较 大 降水 
( 雪 ) ,最 终 使 土壤 水 转变 为 冻 土 ,直至 次 年 接近 融雪 期 ,土壤 水 再 次 增加 且 到 融雪 过 程 结束 。 一 方面 证 明 CMADS+SWAT 模式 
可 有 效 提高 SWAT 水 文 模型 在 我 国 西北 干旱 区 (站 点 稀缺 区 域 ) 的 表现 能 力 ,为 一 方面 理 清 了 精 博 河流 域 相关 地 表 分 量 ( 土壤 
湿度 .蒸发 等 ) 时 空 演 变 规律 。 本 研究 对 我 国 大 气 水 文学 种 发 展 将 起 到 一 定 的 科学 促进 作用 。 
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Abstract: Describing the changing rules of land surface variables by using meteorological and hydrological models has 
always been/@ research hotspot in the fields of atmospheric and hydrological science. However, due to the scarcity of weather 
stations, in West China, traditional weather stations cannot satisfy the requirements of high precision simulation of large scale 
land surface, components. Using the Jing and Bo River Basin in Xinjiang as a research area, the present study used the 
China Meteorological Assimilation Driving Datasets for the SWAT model ( CMADS) to drive the Soil and Water Assessment 
Tool. (SWAT) model, and then completed calibration, verification, and time-space relation analysis of each land surface 
component (such as soil moisture content, snow depth, and snow melt). Our analysis showed that the CMADS dataset can 


drive and calibrate the SWAT mode land complete localization work in Jin and Bo River Basin. The NSE efficiency 
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coefficients of the SWAT model driven by the CMADS dataset were generally controlled between 0.659 and 0.942 at the 
monthly scale, and were also controlled between 0.526 and 0.815 at the daily scale. Furthermore, our analysis indicated 
that the soil moisture content would reach a high level for the first time between March and April each year, which 1s mainly 
caused by snow melting in the high mountains. However, when snow melting finished, due to an increase in precipitation 
and temperature, soil temperature fluctuated until the middle of October, when cold air brought considerable precipitation 
and snow. Finally, soil water was transformed into frozen soil until the snow melting period in the following year. Thereafter, 
soil water would increase again until the end of the snow melting period. On one hand, this study verified that the CMA DS+ 
SWAT mode can enhance the performance ability of the SWAT model in the arid areas of northwestern China, which lacks 
weather stations. On the other hand, the study provided a scientific explanation for time-space changing rules of land surface 
components ( such as soil moisture and evaporation) in the Jing and Bo River Basin. The findings of this research will play a 


certain role in promoting the development of China's meteorological and hydrological sciences. 


Key Words: CMADS; SWAT; Jing and Bo River Basin; surface component 


近年 来 ,全球 范围 内 均 面 临 严 重 的 生态 与 水 危机 ,干旱 区 陆 面 过 程 变化 无 论 对 整个 内 陆 河 流水 循环 亦 或 
是 对 整个 干旱 区 生态 环境 演变 均 具 重要 影响 。 其 中 ,干旱 内 陆 区 水 量 这 k 山 加 速 又 减 与 恶化 也 促使 该 地 区 突 
发 性 水 事件 频繁 发 生 。 干旱 区 绿洲 构造 独特 ,水 与 能 量 循环 具备 其 自身 的 特殊 规律 ,系统 分 析 新 疆 干 旱 区 地 
表 水 资源 相关 分 量 时 空 变化 意义 重大 ,这 将 为 我 国 实现 内 陆 干 旱 区 流域 生态 水 文 可 持续 发 展 提供 重要 科学 依 
据 。 精 博 河 流域 位 于 新 疆 地 区 西部 ,封闭 流域 总 面积 达 25 万 km*。 流 域 在 西北, 南 等 3 个 方向 均 为 高 山 ,中 
间 为 谷地 ,东边 为 例 地 。 地 貌 及 地 形 宏观 上 与 准噶尔 倪 地 连 为 = 条 ,大 体 分 为 山谷 .分 三 大 地 貌 。 流 域 交汇 
处 正好 位 于 阿拉 山口 大 风 下 沉 通 道 , 故 大 风灾 害 是 该 流 虞 最 大 的 气象 灾害 。 此 外 ,由 于 该 流域 生态 结构 极其 
脆弱 ,导致 研究 区 内 风沙 频繁 . 盐 渍 化 及 攻 漠 化 现象 严重 。 近 年 来 ,由 于 气候 变 暖 加 剧 及 人 类 垦荒 活动 不 断 ， 
导致 流域 内 艾 比 湖人 湖水 量 不 断 减少 ,生态 环境 持续 屎 化 。 据 不 完全 统计 ,日 20 世纪 50 ERE 70 年 代 以 
来 , 精 博 河 流域 内 芯 比 湖面 积 萎 缩 了 近 678km“ 哆 比 湖 储 水 量 共 损耗 近 23 亿 m 。 水 资源 的 急剧 下 降 导 致 艾 
比 湖水 体 矿 化 度 急剧 攀升 ,湿地 面积 嗓 然 减 小 ,湖泊 对 气候 的 调控 作用 降低 ,对 流域 内 生态 系统 及 当地 居民 的 
生存 环境 造成 了 极 大 的 破坏 影响 4 并 最 终 导 致 精 博 河 生态 环境 迅速 恶化 ,发 生 了 生物 链 断 层 多样 性 减少 ,元 
省 化 过 程 加 剧 等 现象 。 由 于 精 博 河流 域 对 新 疆 生 态 平 衡 及 社会 经 济 发 展 均 具 有 重大 页 献 ,利用 高 精度 大 气 数 
据 驱 动 成 熟 水 文 模型 的 输 乔 结果 可 为 流域 提供 一 定 的 可 持续 发 展 决 案 , 促 进 精 博 河 流域 地 表 过 程 分 量 ( 和 后 
态 、 水 文 等 ) 合理 发 展 。 

由 于 精 博 河流 域 传 统 和 气象 站 点 稀缺 ,流域 下 热 面 分 异 较 大 ,加 之 气候 变化 及 人 为 活动 的 影响 ,增强 了 该 区 
域 模 式 输出 与 评估 的 不 确定 性 。 众 多 研究 证 明 : 融 合 更 多 观测 数据 的 大 气 同 化 数据 集 可 更 大 程度 提高 模 
式 输出 精度 坚 3 。 当 前 ,国内 外 大 气 再 分 析 种 类 繁多 ,如 : NCEP 的 NCAR/DOE'*". ,欧洲 预报 中 心 ERA 一 15、 
ERA—40.,ERA--Interim 再 分 析 资 料 551 日 本 JRA—25 再 分 析 资 料 ' 引 及 陆 面 模式 驱动 场 Princeton 数据 集 。 由 
于 传统 火气 观测 人 站 并 未 布 满 全 球 , 以 上 再 分 析 数 据 集 对 气象 缺 测 区 气候 场景 分 析 及 应 用 提供 了 重要 数据 支 
BE. 然而 ,对 于 区 域 精细 化 模拟 而 言 ,以 上 再 分 析 数 据 的 准确 性 示 待 提高 。Pal 等 … 利用 区 域 气候 模式 
RegCM3 对 东亚 季风 区 冬 . 夏 两 个 季节 内 降水 量 的 月 变化 规律 在 季节 尺度 上 进行 了 模拟 与 评估 后 发 现 ; 区 域 
模式 RegCM3 对 降水 量 存在 较 大 误差 ,这 种 现象 在 冬季 显得 尤为 凸显 。 赵 天 保 等 利用 NCEP ERA 及 JRA 等 
数据 集 的 地 面 压强 及 温度 要 素 在 中 国 区 域 可 用 性 进行 评估 后 发 现 :3 种 再 分 析 数 据 集 存 在 明显 的 季节 性 和 区 
域 性 差异 中。 施 晓 上 晕 等 “利用 各 类 技术 手段 将 NCAR 驱动 数据 (包括 :空气 温度 风速) 在 中 国 区 域 的 适 
用 性 进行 评估 后 发 现 : 风 场 数据 距 平 误差 与 海拔 表现 为 负 相 关 现 象 。 以 上 人 研究 均 表 明 . 各 类 再 分 析 产 品 虽 可 
较 好 地 反映 出 大 尺度 气象 要 素 空 间 模 态 ,然而 ,由 于 大 部 分 再 分 析 数 据 集 并 未 利用 我 国 全 部 区 域 自动 站 进行 
订正 , 故 不 能 较 好 地 反映 我 国 近 地 表 真实 气象 状况 。 综 上 所 述 : 利 用 更 接近 我 国 真实 气象 场 的 大 气 数据 集 驱 
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动 成 熟 的 水 文 模型 ,可 更 好 的 对 我 国 大 尺度 地 表 分 量 时 空 演变 规律 进行 细致 化 分 析 | 。 

本 人 研究 以 站 点 稀缺 .下 垫 面 分 异 较 大 的 新 疆 博 精 河流 域 作为 典型 验证 及 分 析 区 ,利用 SWAT 模型 中 国 大 
气 同化 驱动 数据 集 (CMADS) ”驱动 在 我 国 应 用 普遍 的 SW AT 水 文 模型 ,以 期 在 为 该 流域 提供 有 效 生态 水 文 
系统 恢复 策略 的 同时 ,验证 CMADS 数据 集 在 中 国 西部 地 区 的 适用 性 。 


1 研究 区 介绍 


精 博 河流 域 地 处 西天 山北 坡 ( 图 1) , 4 F 81946'—83?51'E,44?02'—45?10'N 之 间 。 总 面积 达 11275km?, 
流域 内 降水 多 为 来 目 北 冰 洋 和 大 西洋 的 水 汽 ,呈现 总 体 山区 较 平 原 多 ,西部 较 东 部 多 , 阴 坡 较 阳 坡 多 的 态势 。 
精 河 和 博 尔 塔 拉 地 区 的 冰川 面积 均 分 别 为 96.2km* 和 110.3km*。 其 中 , 精 河 冰川 补给 量 在 年 际 考 约 为 0.96 亿 
m , 博 尔 塔 拉 河 约 为 1.05 亿 m ,分 别 占据 了 其 本 身 河川 总 径流 的 20.6% 和 21.4%。 近 年 来 ,流域 内 点 比 湖 湖 
水 矿 化 度 在 逐年 升 高 ,流域 生态 与 生活 用 水 矛盾 也 随 之 日 益 凸 显 。 流 域 土地 退化 现象 严重 ,其 中 ;区 比 湖 流 域 
湖底 已 有 1500km ”区域 退化 沦 为 盐 漠 , 盐 渍 化 面积 达 71km 。 
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1 博 、 精 河流 域 示 意图 


Fig.1 Schematic diagram of the Jing and Bo River Basin 


Yr A3 , "r8 DORT RS E387 i398 77 E 8 ^p 39 88 DC BAL EA UE 2 DU T 8 2S BEA EOS (3X DORIA, ) 及 
ERARI. HTAA 4 个 气象 站 ,导致 对 该 区 域 还 未 进行 较为 合理 可 信 的 水 文 过 程 模拟 与 分 
析 工 作 。 本 研究 将 基于 面 回 东亚 的 高 分 辨 卒 大 气 驱 动 场 CMADS 数据 集 驱 动 SWAT 模型 ,对 精 博 河水 文 相关 
分 量 进行 细致 模拟 ,以 期 系统 分 析 精 博 河流 域 各 地 表 要 素 时 空 演变 特征 ,同时 为 我 国 高 分 辩 率 气象 数据 驱动 
FAIRER ATT T E SEAE o 


2. 材料 与 方法 


2.1 SWAT 模型 地 表 输 入 数据 
2.1.1 数字 高 程 模型 数据 

SWAT 模型 输入 数据 主要 包括 数字 高 程 模型 (DEM ) 土壤 数据 、 土 地 利用 数据 及 气象 驱动 场 。 其 中 , 数 
字 高 程 模 型 选取 自 SRTM —(90m) DEM( FÉ 2) , 该 数据 来 源 于 CGIAR—CSI SRTM 90m 高 程 数 据 库 ( http:// 
srtm.csi.cgiar.org/ SELECTION/ inputCoord.asp ) 。 
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2 研究 区 高 程 及 SWAT 模型 其 他 输入 数据 
Fig.2 The elevation of the study area and the input data for SWAT model 


2.1.2 ” 精 博 河流 域 土壤 分 布 

土壤 物理 特性 决定 了 SWAT 模型 不 同 水 文 响应 单元 的 产 时 世 流 特征 ,同时 也 是 SWAT 模型 定义 水 文 啊 应 
单元 (HRU ) 的 参考 依据 。 本 研究 选用 的 土壤 输入 数据 内 基于 世界 土壤 数据 库 (HWSD ) 的 中 国土 壤 数 据 集 
(V1.1) (图 3A)'”。 土 地 利用 来 源 于 精 博 河流 域 管理 处 ,土地 利用 数据 又 加 了 我 国 第 二 次 冰川 数据 集 "” Lee 
加 后 的 土地 利用 分 布 情况 可 参见 图 3B。 为 保证 模式 分 辨 率 一 致 ,本 人 研究 将 DEM .土壤 与 土地 利用 等 数据 空 
间 分 状 率 均 统一 为 Tm ,平面 投影 坐标 统 生 设 顽 为 WGS_1984_UTM _Zone_44N。 
2.1.3 ”大气 驱动 数据 

本 研究 使 用 的 大 气 输入 数据 SWAT 模型 中 国 大 气 同 化 驱动 数据 集 ( China Meteorological Assimilation 
Driving Datasets for the SWAT imodel, CMADS) ,该 数据 集 引 入 了 中 国 陆 面 同化 系统 ,保证 了 气象 输入 数据 的 可 
HEE ”。CMADS Xl ARAR B Dose 1 KK 2, 


x1 CMADS 数据 集 信 息 
Table 1 The information of CMADS driving data 


数据 集 CMADS ( V1.0) 
Datasets China Meteorological Assimilation Diving Datasets for the SWAT model ( Version 1.0) 
提供 要 素 Provide elements 日 最 高 / 低 气 温 .日 平均 相对 湿度 .日 均 风 速 .日 累计 降水 量 .日 累计 太阳 辐射 
原始 空间 范围 Original space range 0—65?0'N,60—160?0'E 
本 研究 空间 范围 /Spatialgcope 81°46'—83°51'E  44202'—45?10'N 
Ehi NL E Time scale 2008.1.1—2014.12.31 
原始 分 辨 率 Original resolution 1/3?,1/4?,1/8?,1/16? 
本 人 研究 使 用 空间 分 辨 率 Spatial resolution in this study 1/3? 
站 点 数 Station number 23 个 
CMADS V1.0 版 本 空间 分 辨 率 为 1/3 度 , 时 间 分 辨 率 为 逐日 ,数据 尺度 为 2008 一 2014 年 。SWAT 模型 在 


精 博 河流 域 共 读 取 了 CMADS 数据 集 近 23 个 站 点 要 素 ( 温 湿 、 风 、 降 水 及 辐射 数据 ) 。 其 中 ,CMADS 气象 站 
点 具体 信息 (站 点 名 称 及 经 纬度 ) 如 表 2 所 示 。 
2.1.4 水 文 验 证 数据 

本 研究 将 精 博 河流 域内 水 文 站 逐日 流量 资料 作为 模型 率 定 及 验证 资料 , 台 站 信息 如 表 3 所 示 。 
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图 3 精 博 河流 域 土壤 分 布 及 精 博 河流 域 土地 利用 分 布 
Fig.3 Soil distribution in the Jing and Bo River Basin, and Land use distribution in the Jing and Bo River Basin 
iE :SWGR : 细 芽 草 , Slender Wheatgrass ; PAST ; 1437; , Pasture; HAY ; TF , Hay; WETL : HA IEH , Wetlands-Mixed ; WATR : 7K 3 , Water; AGRL ; H 
他 作物 , Agricultural Land-Generic ; AGRR ; "P SEJEJ , Agricultural Land-Row Crops; URML :居民 区 , Residential-Med/Low Density; UIDU ; 工业 用 
Ji , Industrial ; UTRN :交通 用 地 ,Transportation; FRST ; 林地 , Forest-Mixed; RNGE ; 草场 , Range-Grasses; BARO ; ZR Z1, Bare rock; DESE: 荒漠 ， 
Desert; SALA :盐碱地 ,Saline land; GLSN :冰川 /永久 积 雪 ,Glacier and snow; FRSD :落叶 林 , Forest-Deciduous ; PINE ; £| HER, Pine; RNGB :灌木 
Jk, Range-Brush ; ORCD :果园 , Orchard ; FRSE : 常 绿林 ,Forest-Evergreen 
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表 2 精 博 河流 域 CMADS 站 点 分 布 
Table 2 The station distribution of CMADS in Jing and Bo River Basin 


eet 纬度 dm 经 度 EC 高 程 /m eke 纬度 NC) 经 度 高 程 /m 
of CMADS Latitude Longitude Elevation of CMADS Latitude Longitude Elevation 
136—62 44.99 80.34 2968.00 136—68 44.99 82.34 583.00 
136—65 44.99 81.34 1074.00 135—67 44.65 82.01 504.00 
136—69 44.99 82.68 184.00 135—68 44.65 82.34 296.00 
137—967 45.32 82.01 1863.00 135—64 44.65 81.01 2072.00 
136—66 44.99 81.68 735.00 135—606 44.65 81.68 1432.00 
136—71 44.99 83.34 232.00 134—69 44.32 82.68 936.00 
136—967 44.99 82.01 649.00 134—68 44.32 82.34 1652.00 
136—70 44.99 83.01 206.00 134—67 44.32 82.01 2081.00 
136—72 44.99 83.67 249.00 134—70 44.32 83.01 1898.00 
135—71 44.65 83.34 236.00 134—71 44.32 83.34 3057.00 
135—70 44.65 83.01 206.00 133—70 43.99 83.01 2338.00 

135—060 44.65 82.68 263.00 


CMADS 站 点 名 为 CMADSV1.0 数据 集 在 其 空间 范围 (东亚 区 域 ) 内 唯一 标识 符 


表 3 精 博 河流 域 水 文 站 相关 信息 


Table 3 Hydrological stations data statistics in Jing and Bo River Basin 


水 文 站 纬度 经 度 台 站 海拔 /]m 数据 时 段 /a 
Hydrological station Latitude Longitude Station elevation Data interval 
ili 5& Wenquan 44*59'N 81°02% 1310 2009—2013 
精 河 Jinghe 44?22'N 82?*55'E 620 2009—2013 


2.0 SWAT 水 文 模式 介绍 

SWAT 模型 概 化 的 将 相同 土地 柳 被 利 儿 s 二 二 分 布 及 管理 方式 的 均一 性 单元 视 为 同一 类 水 文 响应 单元 
(HRU) ,模型 认为 每 一 种 HRU 内 所 有 水量 平衡 过 程 具 一 致 性 特征 。SWAT 模型 提供 了 3 种 潜在 蒸发 计算 
方法 , 主要 包括 Penman-Monteith-34; 7", Priestley-Taylor 法 及 Hargreaves 法 ~ 。 在 本 研究 中 , 由 于 
CMADS 数据 集 可 提供 全 部 输入 要 素 (而 传统 气象 站 由 于 无 法 提供 太阳 种 射 数据 仅 能 利用 天 气 发 生 器 模拟 )， 
本 研究 选取 Penman-Monteith 活 帮 为 蒸发 模拟 的 方法 ,以 最 大 可 能 提高 模式 异 拟 精度 52 。 


3 SWAT 模型 方案 设置 


本 研究 区 流域 面积 经 过 ArcGIS 提取 后 为 2482366.4052hm2 , 共 划 分 39 个 子 流 域 及 1648 个 水 文 响应 单 
元 。 由 于 降水 输入 数据 为 逐日 数据 ,地 表 径流 模拟 方法 选取 径流 曲线 数 模型 计算 ,地 表 径 流 将 被 分 别 在 不 同 
水 文 啊 雇 单元 电 进 行 模拟 ,最 终 汇流 到 主 河道 。 最 后 ,选取 基于 连续 方程 的 河 槽 苹 量 法 计算 主 河道 水 量 演进 。 
SWAT 和 异型 通过 Centriod 法 将 流域 内 单 点 气象 数据 空间 离散 插值 到 整个 流域 ,为 减少 空间 离散 插值 之 来 的 误 
25 (JU Be SEA DC) ,本 人 研究 针对 精 博 河流 域 划分 出 多 个 高 程 沉 以 期 精确 识别 不 同 高 程 囊 降水 量 分 布 状 态 。 

综合 CMADS 驱动 场 时 间 尺 度 (2008 一 2014 年 ) 及 径流 观测 数据 的 重合 区 间 ,同时 为 使 模拟 初期 所 有 水 文 
过 程 从 初始 状态 进入 平衡 状态 ,本 研究 将 预 热 期 设 定 为 1 年 ( 即 2008 年 ) ,校准 期 设 为 2009 一 2010 年 ,验证 期 
设 定 为 2011 一 2013 年 。 
3.1 敏感 性 分 析 

本 人 研究 利用 SWAT-CUP 对 SWAT 模式 进行 模式 率 定 工作 。SWAT-CUP 是 EWAGE 研究 所 针对 SWAT Bi 
型 研发 的 一 款 自 动 校准 与 不 确定 性 分 析 的 程序 。 在 本 研究 中 ,重点 对 12 个 与 径流 相关 的 参数 进行 参数 率 
E ,并 得 到 模式 参数 最 终 取 值 ( 表 4) 。 
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3.2 ”模型 校准 

在 模型 率 定期 (2009 一 2010 年 ) 对 模型 参数 进行 校准 工作 。 笔 者 先 利 用 SWAT-CUP 率 定 软 件 对 精 博 河 
流域 两 个 测 站 同时 进行 参数 校准 ,校准 过 程 中 兼顾 参考 流域 年 均 莱 发 量 及 降水 量 情况 下 进行 宏观 调整 ,首先 
考虑 年 雁 发 量 和 径流 量 的 关系 比 ,保证 模拟 结 末 符合 精 博 河流 域 年 总 燕 发 量 、 降 水 量 及 径流 量 接近 实际 真实 
情况 。 分 然后 延展 至 月 ,在 月 数据 校准 准确 后 进行 日 参数 率 定 及 微调 。 分 析 发 现 :SWAT 模型 降水 递减 率 
(PLAPS) 参 数 最 优 取 值 为 44.5mm/km。 温 度 梯度 (TLAPS) 最 优 值 为 -4.3039%C/km( 表 4) ,与 流域 内 多 年 实 
味 均 值 吻合 。 


表 4 SWAT 模型 参数 最 终 取 值 结果 
Table 4 The final results of the model parameters within SWAT model 


变量 名 参数 定义 最 终 优 选 值 
Variable Parameter definition Final optimal value 
CN2.mgt SCS 径流 曲线 值 34 
ALPHA, BF.gw 基 流 a BIF (days) 0.453 
GW_DELAY.gw 地 下 水 补给 延迟 时 间 (days ) 42 
GWQMN.gw 地 下 水 汇 和 人 主 河道 时 浅 层 含水 层 的 水 位 国 值 (mm) 8.5 
GW. REVAP.gw 地 下 水 再 蒸发 系数 0.03 
ESCO.hru 土壤 蒸发 补偿 系数 1.008 
ALPHA_BNK.rte 基 流 退 水 常数 (days) 0.45 
SFTMP 降雪 日 的 平均 空气 温度 (Y ) 4.8 
PLAPS 降水 直 减 率 (mm/km) 44.5 
SMFMN 12 月 21 日 的 融雪 因子 (mm/d-%C) 9.407 
SMFMX 6 H 21 日 的 融雪 因子 (mm/d-%C) 0.1302 
TLAPS 温度 直 减 率 ( /km) -4.3039 


CN2:SCS 径流 曲线 值 , Initial SCS runoff curve number for moisture condition II; ALPHA, BF; 3&jj a 因子 ,Baseflow alpha factor; GW. DELAY ; 3t F 
水 补给 延 退 时 间 , Groundwater delay time; GWQMN : JH F KRC A E ANE EST 26 E 2K JB ZK f B] ME , Threshold depth of water in the shallow aquifer 
required for return flow to occur; GW.. REVAP ; 4E F ZK-HEz& A AR, Groundwater " revap" coefficient; ESCO ; 土壤 蒸发 补偿 系数 , Soil evaporation 
compensation factor; ALPHA, BNK ; 3£ ji; 3H 7K 2$ Zt; Baseflow alpha factor for bank storage; SFTMP ; 降雪 日 的 平均 空气 温度 , Snowfall temperature; 
PLAPS :降水 直 减 率 ,Precipitation lapse rate; $óMFMN; 12 H 21 日 的 融雪 因子 ,Melt factor for snow on December 21; SMFMX:6 H 21 日 的 融雪 因子 ， 
Melt factor for snow on June 21; TLAPS :温度 直 减 率 ,Temperature lapse rate 


3.3 ”模型 评估 

本 研究 选取 (NSE Nash-Sutcliffe Efficiency ) 效率 系数 及 R 决定 性 系数 作为 模型 评估 指标 ,以 上 指标 被 广 
泛 的 应 用 于 评估 模型 的 性 能 站 。 其 中 ,Nash-Sutcliffe 效率 系数 是 一 个 正 态 统计 方程 , 它 反 映 了 观测 值 与 对 
应 模拟 值 的 拟 舍 程度 @ NSE 可 利用 方程 (1) 计算 一 二) ， 


2.50809 
NSE = 1 -一 — (1) 
2. (Qni a Qna) 
WEF ,万 为 径流 值 变量 , 0, M 0. 分 别 代表 观测 值 及 模式 值 。0, 表示 径流 平均 观测 值 。NSE 变化 范围 为 -oo 至 
1,25 NSE 计算 结果 为 1 时 ,认为 观测 值 与 模拟 值 完全 吻合 ; 当 其 值 为 0.5 一 1 之 间 时 ,模式 结果 可 接受 , 当 NSE 
小 于 0 时 ,认为 模式 结果 较 差 。 
决定 性 系数 RR 表示 变量 之 间 相 关 程 度 ( 见 式 2) : 


[3(09,,-02(0,; -0)] 
ccc 一 - (2) 
2 (070. "X Oeo 


i 
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式 中 , Q, , Q. 和 式 (6) 中 的 含义 相同 ,i 变量 表示 第 i 次 的 模拟 或 实测 值 。 

众多 研究 将 R 20.5 H. NSE>0.5 作为 SWAT 模型 模拟 满意 的 准则 ” , 亦 有 研究 认为 NSE>0.4 也 可 作为 
模型 指标 满意 的 标准 ,本 研究 采用 Moriasi 等 评价 标准 [9 。 即 模型 在 率 定 期 的 月 尺度 模拟 结果 NSE 20.65 或 
日 尺度 模拟 结果 NSE 20.5 ,说明 模型 模拟 结果 可 接受 。 


4 结果 分 析 


41 基于 CMADS+SWAT 模式 的 征 流量 模拟 验证 及 分 析 率 定期 
本 研究 利用 CMADS+SWAT 模式 分 别 输出 流域 两 
内 两 个 水 文 控制 站 ( 精 河山 口 控制 站 及 温泉 站 ) 逐 月 径 
流量 ,利用 观测 结果 进行 参数 率 定 与 校准 工作 (图 4 一 
图 7)。 
分 析 发 现 :在 月 尺度 上 ,由 CMADS 驱动 下 的 SWAT 
模型 模拟 结果 在 精 博 河流 域 两 个 控制 站 均 到 达 了 满意 


验证 期 A 95PPU 
A 观测 值 
A 模拟 值 


EN 


o2 


径流 量 Flow out/(m?/s) 


10 20 30 40 50 60 
的 效果 ( 见 表 5)。 由 CMADS 驱动 下 的 模拟 结果 在 精 河 时 间 Time/months 
山口 站 的 NSE 和 R? 均 取得 了 令 人 满意 的 表现 (NSE = 4 CMADS 数据 集 驱动 SA 了 模型 的 精 河山 口 控制 站 逐 月 模 
0.939, R^ 20.942) , MER (2009—2013), dB ,95PPU 代表 预测 不 确定 性 内 监测 数 


在 模型 验证 期 ,虽然 NSE 效率 系数 与 确定 性 系数 EUER 

REER ss HH p ME 也 得 到 了 可 接受 的 效果 ( NSE - 0. Fig.4 The ae and MEN nn d NN 
by SWAT. model driving by CMADS in Jinghe control station, 

904 ,R" 3 0.934) S 相对 精 河 山 H 控制 站 来 次 2 LES 站 在 where the percentage of measured data bracketed by the 95% 

率 定 期 和 验证 期 的 模拟 精度 均 较 前 者 KR, 笔者 认为 , prediction uncertainty ( 95PPU) 

ie Ef eg EL KAKI DE e P AZ ARS] ELE, , EH 

T SWAT 模型 度 日 因子 仅 对 融雪 进行 了 细致 化 考虑 ,导致 流域 在 温 录 站 的 模拟 精度 要 低 于 精 河 控制 站 (冰川 

补给 率 较 温 果 控制 站 小 ) 所 在 子 流域 。 


加 率 定期 (月 ) TM 验证 期 (月 ) 
a b 

45 y= 0.9976x - 0.0007 40 + y= 1.094x - 0.3648 

40 R? = 0.9423 35 上 R? = 0.9342 


径流 量 观 测 值 


Runoff observed/(m?/s) 


0 
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径流 模拟 值 
Runoff simulated/(m?/s) 


5 精 河山 口 控制 站 月 观测 与 模式 模拟 结果 确定 性 系数 


Fig.5 The certainty factor of the relationship between monthly observation and model simulation data in Jinghe control staion 


当 完 成 月 尺度 率 定 和 验证 后 ,本 研究 将 月 尺度 最 佳 参数 值 带 入 SWAT 模型 进行 日 微调 校准 后 发 现 . 
CMADS+SWAT 模式 径流 模拟 结果 在 两 个 水 文 控制 子 流域 与 日 观测 结果 呈现 很 好 的 一 致 性 。( 兄 网 8 一 图 
11、 表 5)。 

在 模型 率 定期 ,由 CMADS 驱动 SWAT 模式 在 精 河 站 的 日 拟 合 结 果 ( NSE 20.801, R^ 2 0.815) 5i if 5 vi 0 
合 结果 (NSE= 0.796, 尺 =0.791) 也 均 较 好 。 在 验证 期 ,虽然 SWAT 模式 在 两 个 站 点 模拟 结果 均 可 接受 ,而 模 
式 在 精 河山 口 站 模拟 效果 (NSE=0.851 , R^ 20.796) $i 5 Ta tll yi ( NSE- 0.526, R^ 20.592) 表现 更 佳 。 
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4.0 基于 CMADS+SWAT 模式 的 土壤 湿度 及 相关 变量 分 


将 SWAT 模型 进行 精 博 河流 域 本 地 化 后 ,为 从 时 空 
角度 分 析 CMADS+SWAT 模式 在 土壤 湿度 及 融雪 变量 
时 空 演 变 模 拟 能 力 , 同 时 定量 分 析 各 分 量 之 间 响 应 关 
系 , 本 研究 还 以 精 河 子 流域 为 主要 分 析 对 象 ,对 多 个 要 
素 啊 应 关系 进行 了 提取 。 

4.2.4. 精 博 河流 域 融 雪 过 程 与 土壤 湿度 啊 应 分 析 

为 俩 究 精 博 河流 域 融 雪 过 程 对 土壤 湿度 的 影响 ,本 
部 分 提取 精 博 河 整 个 流域 在 2012 年 4 月 8 日 的 土壤 湿 
度 及 其 对 应 融雪 量 空间 变化 图 (图 12) 。 

图 12A 为 当日 土壤 温度 分 布 情况 ,图 12B 为 对 应 
时 期 融雪 空间 分 布 情况 。 为 定量 研究 流域 内 各 地 表 分 
量 间 关系 ,研究 提取 了 精 河 控制 站 子 流域 内 各 地 表 分 量 
(包括 土壤 温度、 潜在 苹 散 发 .降水 量 和 融雪 量 ) 在 2012 
年 4 月 8 日 的 时 间 序 列 分 布 (图 13)。 
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6 CMADS 数据 集 驱 动 SWAT 模型 的 温泉 控制 站 逐 月 模拟 结 
果 (2009 一 2013)， 其 中 ,95PPU 代表 预测 不 确定 性 内 监测 数据 占 
据 的 百分比 

Fig.6 The monthly simulation results simulated by SWAT model 
driving by CMADS in Wenquan control station, where the 
percentage of measured data bracketed by the 95% prediction 
uncertainty ( 95PPU) 
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Fig.7 The certainty factor of the relationship between monthly observation and model simulation data in Wenquan control station 


通过 对 2012 4E 44H 8 日 精 情 河流 域内 土壤 湿度 分 
布 分 析 发 现 : 流域 整体 填 壤 温度 在 该 日 处 于 侦 湿 状态 
(图 14) 。 其 中 , 温 果 站 控制 子 流域 土壤 湿度 为 18mm 
至 20mm 左右 , 而 精 河 控制 站 及 其 附近 子 流域 土壤 湿度 
达 近 57mm, 遂 过 分 析 流 域内 对 应 融雪 空间 分 布 情况 
(图 14) 后 友 现 ,2012 年 4 月 8 日 , 精 博 河流 域内 出 现 整 
体 较 头 融 必 现 绷 , 而 发 生 融 雪 的 区 域 恰 好 为 西天 山北 坡 
地 区 ,该 区 域 年 内 3 一 4 月 雪 盖 覆盖 丰厚 。 为 定量 分 析 
精 博 河 流域 内 融雪 及 土壤 湿度 量 级 及 其 直接 响应 关系 ， 
本 人 研 究 对 精 河 水 文 站 控制 子 流域 中 土壤 湿度 与 融雪 量 
等 分 量 进 行 了 (以 精 河 山口 控制 站 为 例 ) 时序 量 级 分 析 
人 研究 。 

图 13 给 出 了 精 河 控制 子 流域 内 各 地 表 或 近 地 表 分 
量 ( 洪 在 蒸 散发 ,土壤 湿度 ,降水 量 及 融雪 量 等 ) 在 精 河 
山口 控制 子 流 域 时 序 分 布 状 态 。 分 析 发 现 : 精 博 河流 域 
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8  CMADS 数据 集 驱动 SWAT 模型 的 精 河山 口 控制 站 逐日 模 
拟 结果 (2009 一 2013) ,其 中 ,9SPPU 代表 预测 不 确定 性 内 监测 数 
据 占据 的 百分比 

Fig.8 The daily simulation results ( 2009—2013) simulated by 
SWAT model driving by CMADS in Jinghe control station, where 
the percentage of measured data bracketed by the 95% prediction 
uncertainty ( 95PPU) 
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Fig.9 The certainty factor of the relationship between daily observation and model simulation data in Jinghe.control staion 


H 2012 年 3 月 中 旬 开 始 出 现 融雪 现象 ,并 于 2012 年 4 
月 8 日 出 现 较 大 融雪 现象 。 其 中 , 精 河 站 控制 子 流域 融 
雪 量 达 18mm/d。 当 日 精 河 控制 站 子 流域 土壤 湿度 也 
达到 其 年 内 较 高 水 平 ( 约 56mm) 。 分 析 表 明 : 融 雪 期 
内 , 精 博 河流 域 土壤 湿度 增加 的 最 主要 贡献 来 自 于 融 
雪 , 仅 极 少 量 来 自 降 水 。 此 外 ,土壤 湿度 与 潜在 蒸发 数 
值 上 呈 反 比 关系 。 
4.2.2，” 精 博 河流 域 降水 与 土壤 湿度 响应 分 析 

前 文 提 到 精 河 站 控制 子 流 域 融雪 过 程 对 流域 内 填 
培 湿 度 具 有 较 大 影响 。 笔 者 认 为 降水 量 对 土壤 湿度 页 
献 亦 不 容 忽视 。 笔 者 在 此 重点 分 析 了 夏季 与 秋 来 降水 
对 土壤 湿度 的 影响 。 由 于 精 博 河流 域 土 壤 湿 度 在 融雪 
后 期 至 年 内 9 月 会 出 现 较 长 时 期 波动 ,同时 10 明 中 后 
期 将 发 生 陡 然 上 升 现象 并 保持 恒定 等 现象 ) 笔 者 针对 以 
上 两 阶段 降水 量 与 土壤 湿度 响应 关系 进行 时 空 提 取 分 


AM 率 定期 (日 ) "e 


45 y — 0.8565x + 1.5128 
40 


35 
30 


径流 量 观测 值 
Runoff observed/(m?/s) 


5 15 25 35 45 


m 率 定期 EMT. A 95PPU 
80 | 入 观测 值 
70 | 入 模拟 值 
60 | 


50 


径流 量 Flow out/(m?/s) 


10 &, AMA as. i om ii iha 


i ma di 


200 400 600 800 1000 1200 1400 1600 1800 
日 序 Day 


10 CMADS 数据 集 驱 动 SWAT 模型 的 温泉 控制 站 逐日 模拟 
结果 (2009 一 2013) ， 其 中 ,9SPPU 代表 预测 不 确定 性 内 监测 数据 
占据 的 百分比 

Fig.10 The daily simulation results ( 2009 一 2013 ) simulated by 
SWAT model driving by CMADS in Wenquan control station, 
Where the percentage of measured data bracketed by the 95% 
prediction uncertainty ( 95PPU) 
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Fig.11 The certainty factor of the relationship between daily observation and model simulation data in Wenquan control station 


析 ( 图 14). 2g Scu E LU LET T5 ri] ns E PRESE 2E LEER P RARE Tt CH : UR dC SS DR) 演变 情况 ,研究 还 


提取 了 精 河 山口 控制 站 冻 土 和 雪 深 逐日 台 站 观测 资料 ( 
图 14 给 出 了 精 博 河流 域 土壤 湿度 与 降水 量 在 日 内 
况 。 其 中 ,图 14 为 2010 年 6 月 22 日 精 博 河 流域 土壤 


2010 至 2011 4E) (图 16) 作 为 验证 分 析 。 
(2010 年 6 月 22 日 和 2011 年 10 月 21 日 ) 空 间 分 布 情 
湿度 分 布 情况 ,图 14 为 对 应 流域 内 降水 量 分 布 情况 。 


分 析 发 现 :当日 精 河山 口 控制 站 子 流 域 土壤 湿度 达 24.8mm, 精 河山 口 控制 站 子 流 域 的 降水 量 达 30.1mm。 分 
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12 ” 精 博 河 流域 土壤 湿度 与 融雪 量 空间 相关 分 析 


Fig.12 The correlation analysis of soil moisture and snow melting in Jinghe River Basin 
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13， 精 博 河流 域 各 类 分 量 时 间 序列 分 析 ( 以 精 河山 口 控制 站 为 例 ) 


Fig.13 The correlation analysis in Jinghe River Basin (taking the control station of Jinghe as an example) 
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14 精 博 河流 域 土壤 湿度 与 降水 量 空间 相关 分 析 


Fig.14 The correlation analysis of soil moisture and precipitation in Jinghe River Basin 


析 表 明 :降水 量 在 日 内 大 幅度 增加 导致 该 子 流域 土壤 湿度 在 短 时 期 快速 攀升 。 类 似 情 况 同 时 出 现 于 2011 4E, 
图 14 给 出 了 精 博 河 流域 2011 年 10 月 21 日 土壤 湿度 的 空间 分 布 状态 ,图 14 为 对 应 时 期 降水 量 在 整个 流域 
的 分 布 状态 。 该 时 有 段 精 河 控 制 站 及 温 果 控制 站 所 在 的 子 流域 土壤 湿度 较 大 ,其 中 温 果 站 控制 站 所 在 子 流域 土 
HERBES 18mm 左右 , 精 河山 口 控 制 站 所 在 子 流域 土壤 湿度 达到 近 22mm。 分 析 还 发 现 , 在 流域 内 土壤 湿度 达 
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到 该 时 期 较 高 水 平 的 同时 ,降水 过 程 在 精 博 河 整 个 子 流 域 均 有 发 生 。 为 准确 的 验证 土壤 湿度 与 降水 量 在 量 级 
上 的 响应 关系 ,本 人 研究 提取 了 精 河 山口 控制 站 (2010 年 至 2012 年 ) 其 他 地 表 分 量 ( 湾 在 蒸 散发 土壤 湿度 LE 
水 量 及 融雪 量 ) 对 应 时 序 关系 曲线 (图 15) 。 


表 5 CMDAS+SWAT 模式 径流 输出 的 月 .日 拟 合 结果 评估 
Table 5 The assessment under monthly and daily scale simulated by CMADS+SWAT mode 


时 间 尺 度 精 河 山口 控制 站 Jinghe control station 温 果 控制 站 Wenquan control station 

Time scale 率 定期 NS 率 定 期 验证 期 NS 验证 期 率 定 期 NS 率 定 期 验证 期 NS 验证 期 
Calibration NSE Calibration R? Validation NSE Validation R? Calibration NSE Calibration R? Validation NSE Validation R? 

H Month 0.939 0.942 0.904 0.934 0.917 0.914 0.659 0.844 

H Day 0.801 0.815 0.802 0.851 0.796 0.791 0.526 0:592 


NSE : 纳什 效率 系数 ，Nash-Sutcliffe Efficiency 


分 析 发 现 ,在 2010 年 6 月 22 日 ,流域 内 精 河 山口 控制 站 子 流 域 降 水 量 接近 30mm , 流域 内 土壤 湿度 亦 擎 
JF 28 24.8mm( KI 15) , 2011 年 10 H 21 日, 精 河山 口 控制 站 所 在 子 流域 内 降水 量 达 近 和 mm ,于 此 同时 , 子 流 
域内 土壤 湿度 值 也 快速 攀升 至 22mm 左右 ,土壤 湿度 在 此 后 保持 恒定 状态 (图 15) 。j 秋 未 时 期 , 精 河 山口 控制 
子 流域 土壤 湿度 快速 上 升 现 象 充分 说 明 ,年 末 降 水 对 土壤 湿度 变化 贡献 较 类 。 
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15 ” 精 博 河流 域 土壤 湿度 与 降水 量 时间 序 列 分 析 ( 以 精 河 山口 控制 站 为 例 ) 
Fig.15 The correlation analysis of Soil moisture and precipitation in Jinghe River Basin (taking the control station of Jinghe as an 


example) 


以 年 内 循环 角度 对 土壤 湿度 变化 进行 分 析 后 发 现 :人 研究 区 内 每 年 3 至 4 月 为 融雪 期 ,流域 内 雪 深 经 历 了 
先 上 升 后 下 降 的 过 程 , 并 在 融雪 期 结束 时 雪 深 接近 0m, 这 与 本 研究 在 精 河 站 提取 的 雪 深 观测 数据 一 致 (图 
16) 。 融 雪 期 后 期 气候 转 暖 ,此 时 多 十 期 来 临 ,造成 研究 区 内 土壤 泥 度 快速 葡 升 ,由 于 空气 温度 在 该 时 段 较 
高 , 冻 土 发 生 融化 现象 (图 16) 。 共 发 同时 增加 并 加 速 土壤 水 减少 ,导致 该 时 期 土壤 湿度 呈现 波动 变化 趋势 。 
10 月 后 期 ,流域 由 于 冷 湿 气流 过 境 出 现 大 范围 降水 (降雪 ) ,此 时 土壤 含水 量 将 大 幅 增 加 ,伴随 冷 空 气 来 临 而 
促使 蒸发 量 大 幅度 下 降 。 此 时 ,土壤 液态 含水 量 冻 结 ( 每 年 11 月 至 12 月 ) 并 产生 冻 土 。 季 闻 性 冻 土 将 在 来 
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16 2010 年 .2011 年 精 河 控制 站 雪 深 、 冻 土 观测 值 
Fig.16 The snow depth and frozen soil depth observation in Jinghe control station ( Year 2010-and 2011) 


年 1 一 2 月 到 达 峰 值 直 至 次 年 融雪 期 再 次 陡然 下 降 , 冻 土 再 次 融化 ,A 此 时 ,一 部 份 冻 主将 转换 变 成 土壤 液态 含 
水 量 。 

综 上 所 述 ,流域 内 土壤 湿度 在 每 年 3 一 4 月 份 到 达 其 第 一 次 较 电 水 平 ,这 主要 是 由 于 流域 内 积 雪 、 冻 土 融 
化 而 造成 的 。 融 雪 期 结束 后 ,由 于 降水 量 增多 ,空气 温度 上 天 导致 主 壤 湿度 呈现 上 下 波动 态势 ,直至 10 月 中 
旬 冷 空气 过 境 产 生 较 大 降水 ( 雪 ) 并 最 终 使 土壤 水 转换 为 还 塘 , 至 次 年 融雪 期 ,土壤 液态 水 再 次 增加 直至 融雪 
期 结束 。 


s ”讨论 与 结论 


本 文 利用 CMADS 驱动 SWAT 模型 对 精 精 河 法 域内 两 个 水 文 控制 站 进行 了 径流 精度 模拟 验证 工作 ,利用 
本 地 化 的 SWAT 模型 分 析 了 精 河 控制 站 所 在 子 流域 各 类 地 表 分 量 (包括 土壤 湿度 \ 洪 在 区 散发 .降水 量 和 融 
雪 量 ) 时 空 变化 规律 ,加 深 了 对 该 巴 流 域内 名 分 量 要 系 的 空间 演变 规律 ,解释 了 人 研究 区 内 土壤 湿度 在 不 同 季 
节 的 变化 机 制 。 本 人 研究 后 期 将 继续 考虑 将 更 多 观测 数据 (如 卫星 遥感 数据 ) 同 化 到 SW AT 模型 输出 数据 中 ， 
以 提高 SWAT 模式 的 结果 可 年 性 蛙 在 现 有 阶段 ,利用 CMADS 数据 集 驱 动 SW AT 模型 对 研究 高 分 辨 率 气 象 数 
据 驱 动 复杂 水 文 模型 的 研究 意义 较 大 ,预计 该 数据 集 将 为 我 国 的 水 文生 态 细致 化 模拟 研究 提供 数据 文 掉 。 
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